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1. INTRODUCTION 

Despite the fact that Josephson Junction Arrays (JAA) have been actively stud- 
ied for decades, they continue to contribute to the variety of intriguing and peculiar 
phenomena (both fundamental and important for potential applications) providing at 
the same time a useful tool for testing new theoretical ideas. To give just a few recent 
examples, it is sufficient to mention paramagnetic Meissner effect (PME)[l-5] as well as 
the recently introduced thermophase [6,7] and piezophase [8] effects suggesting, respec- 
tively, a direct influence of a thermal gradient and an applied stress on phase difference 
between the adjacent grains. At the same time, an artificially prepared islands of su- 
perconducting grains [4,5,9,10], well- described by the various models of JJAs, proved 
useful in studying the charging effects in these systems, ranging from Coulomb block- 
ade of Cooper pair tunneling and Bloch oscillations [11,12] to propagation of quantum 
ballistic vortices [13]. 

The present paper reviews some of the recently suggested novel effects [14-16] which 
should manifest themselves either in weak- links-bearing superconductors or in an ar- 
tificially prepared JJAs. In Section 2 we consider the appearance of an electric-field 
induced magnetization (analog of the so-called magnetoelectric effect) in a model gran- 
ular superconductor. The dual effect, that is an appearance of a magnetic field induced 
electric polarization is studied in Section 3. Finally, in Section 4 we discuss a possibility 
of two other interesting effects which are expected to occur in a granular material under 
mechanical loading. Specifically, we predict the existence of stress induced paramag- 
netic moment in zero applied magnetic field (Josephson piezomagnetism) and study its 
influence on a low-field magnetization (leading to a mechanically induced PME). 



2. MAGNETOELECTRIC EFFECT 

To account for the unusual behavior of the critical current under the influence of 
a strong enough applied electric field observed in high-T c superconductos (HTCS) [17- 
20], a possibility of the superconducting analog of the so-called magnetoelectric effect 
(MEE) in JJAs has been recently suggested [15,21]. This effect is similar to (but phys- 
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ically different from) the MEE seen in magnetic ferroelectrics (like antiferromagnetic 
BiFeOz) [22]. It has also been discussed in the context of some exotic nonmagnetic 
normal metal conductors, with the symmetry of a mirror isomer [23], and pyroelectric 
superconductors (where the supercurrent passing through a metal of a polar symmetry 
is assumed to be accompanied by spin polarization of the carriers [24]). The effect 
discussed here entails the electric field generation of a magnetic moment (paramag- 
netic in small to diamagnetic for large electric fields) due to superconducting currents 
that circulate between grains. The basic physical reason for the appearance of this 
effect is that the applied electric field induces a magnetization M, that changes its sign 
(M — > — M) under a time-parity transformation. In the case of weakly-coupled su- 
perconducting grains, the phenomenological reason for the time-parity violation comes 
from the fact that the total free energy F for the material exhibiting both magnetic 
and electric properties contains the term oc mn f mn (E, H), with the coefficients a mn ^ 
(here {m,n} = x,y,z). In the standard MEE, the function f mn = E m H n , which leads 
to the corresponding linear effect, which is either an electric- field induced magnetiza- 
tion M m (E) = dF/dH m = a mn E n (in zero magnetic field) or a magnetic- field induced 

— * 

polarization P m (H) = dF/dE m = a mn H n (in zero electric field). Since H (and M) 
changes its sign under time-parity transformation while E (and P) remains unchanged, 
an electric field induced magnetization will break time-parity symmetry even in zero 
magnetic field applied. As we show below, in our case the symmetry breaking term 
in F, represented by a nonzero coefficient a mn , has a more general nonlinear form for 
fmn- We recall that the standard linear MEE can appear when an external electric field 
E interacts with an inner magnetic field Iidm of the Dzyaloshinskii-Moriya (DM) type 
[25]. The DM interaction leads to a term, apart from the standard isotropic term in the 
Heisenberg Hamiltonian, of the form Hdm = Dij-(SiASj), where Si is a Heisenberg 
spin and the constant vector D it j arises from the spin-orbit coupling. An analogous 
situation occurs in our case, as we describe below. To see how we can get a nonzero a, 
or equivalently a DM type interaction in a granular superconductor, we model a HTCS 
ceramic sample by a random three-dimensional (3D) overdamped Josephson junction 
array. This model has proven to be useful in describing the metastable magnetic prop- 
erties of HTCS [26,27]. In thermodynamic equilibrium, this model has a Boltzmann 
factor with a random 3D-XY model Hamiltonian. Specifically, the general form of the 
Hamiltonian (describing both DC and AC effects) reads 



Here {i} = fj is a 3D lattice vector; J(r it j) is the Josephson coupling energy, with 
= ?i — fj the separation between the grains; the gauge invariant phase difference is 
defined as 



where 4>ij{0) = 4>i — <f)j with 0, being the phase of the superconducting order parameter; 
A it j(t) is (time-dependent, in general) frustration parameter, defined as 



w(0 = E J ( r «)[i-«»^(t)] 




&j(*) = &j(0)-Aj(*), 



(2) 




(3) 
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with A(f, t) the (space-time dependent) electromagnetic vector potential which involves 
both external fields and the electric and magnetic possible self- field effects (see below); 
$o — h/2e is the quantum of flux, with h Planck's constant, and e the electronic charge. 
Expanding the cosine term, and using trigonometric identities we can explicitly rewrite 
the above Hamiltonian as H = — J2i,j J[cx>s(Aij)Si ■ Sj — sm( Aij)k ■ Si A Sj], where the 
two-component XY spin vector is defined as Si = (cos0j, sin0j), and A; is a unit vector 
along the z-axis [28]. We see that the second term in this Hamiltonian has the same 
form as in the DM contribution, and thus we can surmise that the time parity will be 
broken by applying an external field (electric or magnetic) to the granular system. To 
bring to the fore this possibility, we show below in a simple but yet nontrivial model 
that this is indeed the case. 

There are different types of A i:j randomness that can be considered [27]. For simplic- 
ity, in the present paper, we consider a long-range interaction between grains (assuming 
J(rjj) = J) and model the true short-range behavior of a ceramics sample through the 
randomness in the position of the superconducting grains in the array (using the expo- 
nential distribution law f r (fij), see below). Here we restrict our consideration to the 
case of an external electric field only but it can be shown that the scenario suggested 
will also carry through when applying an external magnetic field (which will induce 
another time-parity breaking phenomenon in the granular material; namely, magnetic 
field induced electric polarizability, see the next Section). Besides, in what follows we 
also ignore the role of Coulomb interaction effects assuming that the grain's charging 
energy E c <^ J (where E c = e 2 /2C, with C the capacitance of the junction). 

In the case of a granular material, we show here that the corresponding time-parity 
breaking DM internal field can be related to the electric field induced magnetic moment 
produced by the circulating Josephson currents between the grains. As is known [29,30], 

— * 

a constant electric field E applied to a single Josephson junction (JJ) causes a time 
evolution of the phase difference. In this particular case Eq. (2) reads 

</>ij(t) = <M0) + ^E-f id t (4) 

The resulting AC superconducting current in the junction is 

1p T 

I!j(t) = — sin 0^.(0 (5) 

If, in addition to the external electric field E, the network of superconducting grains 
is under the influence of an applied magnetic field H, the frustration parameter Aij(t) 
in Eq.(3) takes the following form 



7T — * — * 27T — 

A iAt) = ~r( HA R i,j) ■ nj - -f-E ■ f^t (6) 

^0 ^0 



Here, R it j = (r*j + r*,)/2, and we have used the conventional relationship between the 

— * — » — » — » 

vector potential A and a constant magnetic field H = rot A (with dH/dt = 0), as 

— * — * — * 

well as a homogeneous electric field E = —dA/dt (with rotE = 0). In the type 
II HTCS the magnetic self-field effects for the array as a whole are expected to be 
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negligible [31]. The grains themselves are in fact larger than the London penetration 
depth and we must then have that the corresponding Josephson penetration length 
must be much larger than the grain size (since the self-induced magnetic fields can in 
principle be quite pronounced for large-size junctions even in zero applied magnetic 
fields [14]). Specifically, this is justified for short junctions with the size d <C Xj, where 

A j = ^J^o/4:iTfioj c \L is the Josephson penetration length with \l being the grain 
London penetration depth and j c its Josephson critical current density. In particular, 
since in HTCS \l — 150nm, the above condition will be fulfilled for d ~ l//m and 
j c ~ 10 4 A/m 2 which are the typical parameters for HTCS ceramics [29]. Likewise, to 
ensure the uniformity of the applied electric field, we also assume that d <C A#, where 
\ E is an effective electric field penetration depth [20,21]. 

When the AC supercurrent Ifjit) (defined by Eqs.(2), (5) and (6)) circulates around 
a set of grains, that form a random area plaquette, it induces a random AC magnetic 
moment j2 s (t) of the Josephson network [26] 

= E / y(*)(*jAi?y) (7) 

Notice that in the MEE-like effect discussed here for a granular superconductor, the 
electric-field induced magnetic moment in the system is still present in zero applied 
magnetic field due to the phase coherent currents between the weakly-coupled super- 
conducting grains. 

To consider the essence of the superconducting analog of MEE, we assume that 
in a zero electric field the phase difference between the adjacent grains (f>ij(0) = 
which corresponds to a fully coherent state of the array. In this particular case, the 
electric-field induced averaged magnetization reads 

^ T OO 

M S {E) = JJfy = ~JdtJ df itj dR id f(f itj , Rij)i2 s (t), (8) 



where r is the electronic relaxation scattering time, and / is the joint probability 
distribution function (see below). 

To obtain an explicit expression for the electric-field dependent magnetization, we 
consider a site positional disorder that allows for small random radial displacements. 
Namely, the sites in a 3D cubic lattice are assumed to move from their equilibrium 
positions according to the normalized (separable) distribution function 

f(ri,A,j) = ftfj)fR(Rij) (9) 

It can be shown that the main qualitative results presented here do not depend on 
the particular choice of the probability distribution function. For simplicity here 
we assume an exponential distribution law for the distance between grains, f r {r) = 
f(x 1 )f(x 2 )f(x 3 ) with fr(xj) = (l/d)e~ Xj / d , and a short range distribution for the de- 
pendence of the center-of-mass probability fn(R) (around some constant value D). The 
specific form of the latter distribution will not affect the qualitative nature of the final 
result. (Notice that in fact the former distribution function f r (f) reflects a short-range 
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character of the Josephson coupling in granular superconductor. Indeed, according to 
the conventional picture [32] the Josephson coupling J(f%) can be assumed to vary 
exponentially with the distance between neighboring grains, i.e., Jifij) = Je~ K ' Tij . 
For isotropic arrangement of identical grains, with spacing d between the centers of 
adjacent grains, we have k = |, |) and thus d is of the order of an average grain 

size.) Taking the applied electric field along the x-axis, E = (E x , 0, 0), we get finally 

M Z (E X ) = ^1 - H Z (E X ), (10) 
for the induced transverse magnetization (along the x 3 = z-axis), where 



and 



^(^)= M o(^)io g ji+(^y, (12) 



stand for the electric-field induced magnetic induction B Z (E X ) and magnetic field 
H Z (E X ), respectively. The induced Josephson current I(E X ) is simply given by Am- 
pere's law I(E X ) = H z (E x )d. In these equations, M = 2neJNdD/h, with N the total 
number of grains and E = H/2der. Eq.(10) is the main result of this Section, which we 
proceed to analyze below. As is seen from Eq.(10), the behavior of the magnetization 
in the applied electric field is determined by the competition between the two contribu- 
tions, the magnetic induction B Z (E X ) and the current induced magnetic field H Z (E X ) 
(or the corresponding Josephson current I(E X )). Namely, below a critical (threshold) 
field E c 1.94 E where B Z (E X ) > fi H z (E x ), a paramagnetic phase of the MEE is 
due to the modification of the magnetic induction in the applied electric field. On the 
other hand above the threshold field (when E x becomes larger than E c ) the Joseph- 
son current I(E X ) induced contribution starts to prevail, leading to the appearance of 
the diamagnetic signal (seen as a small negative part of the induced magnetization in 
Fig.l). Such an electric field induced paramagnetic-to-diamagnetic transition has been 
actually observed [17-19] for behavior of the critical current in ceramic HTCS and it 
was attributed [20] to a proximity-mediated enhancement of the superconductivity in 
a granular material in a strong enough electric field. Explicitly, the critical current 
in a YBCO sample was found to reach a maximum at E = 4 x 10 7 V/m. To relate 
this experimental value with the model parameters, first of all, we need to estimate an 
order of magnitude of the relaxation time r in a zero applied electric field. This will 
provide an upper limit for the relevant r-distribution in our system. It is reasonable to 
connect zero-field r = r(0) with the Josephson tunneling time [29] tj = (R /R n )(h/ J) 
(where R = h/Ae 2 , and R n is the normal state resistance between grains). Typically, 
for HTCS ceramics J/k B ~ 90K and R n /R ^ 10~ 3 , so that tj ~ 10~ 10 s. At the 
same time, at high enough electric fields where the MEE becomes strongly nonlinear, 
we can expect quite a tangible decrease of the relaxation time. Indeed, for an average 
grain size d ps 1/im, the characteristic field E = 4 x 10 7 V/m (which corresponds to 
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Figure 1: The induced magnetization m(e) = M z /M as a function of the normalized 
applied electric field e = E x /E , according to Eq.(10). The inset shows a blow up of 
the paramagnetic region of m(e) (exhibiting its linear superconducting magnetoelectric 
effect regime). 



the region where a prominent enhancement of the critical current was observed [17- 
19]) introduces a substantially shorter relaxation time t(Eq) = H/2deEo ~ 10~ 16 s, in 
agreement with observations. 

To estimate the relative magnitude of the superconducting analog of the MEE 
predicted here, we can compare it with the normal (Ohmic) contribution to the mag- 
netization 

M n {E) = jl n {t) = - J dt J dfijdRijj '(f id , Rij)fl n , (13) 
o o 

where 

A. = E^«a4') (14) 

Here = Vij/ ' R n is the normal current component due to the applied electric field E, 
with Vij = E-fij being the induced voltage, and R n the normal state resistance between 
grains. As a result, the normal state contribution (for E along the x-axis) reads M n = 
a n E x: with a n = ird 2 DN/R n . Similarly, according to Eq.(10), the low field contribution 
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to the superconducting MEE gives M s ~ a s E x with a s = 2ire 2 J Nr(0)d 2 D/h 2 . Thus, 
at low enough applied fields (when E x <C E ) 

a s r(0) , , 

— ^ , (15) 

where Tj = (R / R n ){%/ J) with R = h/4e 2 . According to our previous discussion 
on the relevant relaxation-time distribution spectrum in our model system, we may 
conclude that r(0) < tj. So, we arrive at the following ratio between the coefficients 
of the superconducting to normal MEEs, namely a s /a n < 1. The above estimate of 
the weak-links induced MEE (along with its rather specific field dependence, see Fig.l) 
suggests quite an optimistic possibility to observe the predicted effect experimentally in 
HTCS ceramics or in a specially prepared system of arrays of superconducting grains. 

We note that in the present analysis we have not explicitly considered the polar- 
ization effects (due to the interaction between the applied electric field and the grain's 
charges) which may become important at high enough fields (or for small enough 
grains), leading to more subtle phenomena (like Coulomb blockade and reentrant-like 
behavior) that will demand the inclusion of charging energy effects in the analysis. 



3. MAGNETIC FIELD INDUCED CHARGING EFFECTS 

This Section addresses a related phenomenon which is actually dual to the above- 
discussed analog of magnetoelectric effect. Specifically, we analyze a possible appear- 
ance of a non-zero electric polarization and the related change of the charge balance in 
the system of weakly-coupled superconducting junctions (modelled by the random 3D 
JJAs) under the influence of an applied magnetic field. Since the field-induced effects 
considered in this Section are expected to manifest themselves in high enough applied 
magnetic fields (with a nearly homogeneous distribution of magnetic flux along the 
junctions) and as long as the Josephson penetration length Aj exceeds the characteris- 
tic size of the Josephson network d (which is related to the projected junction area S, 
where the field penetration actually occurs, as follows S = ird 2 ) : the Josephson current- 
induced "self-field" effects (which are important for a large-size junctions and/or small 
applied magnetic fields) may be safely neglected [33]. Besides, it is known [34] that 
in discrete JJAs pinning (by a single junction) actually concurs with the "self-field" 
effects. Specifically, it was found [34] that the ratio d/\j is related to the dimensionless 
pinning strength parameter (3 as d/\j = ^Jfi suggesting that a weak pinning regime 
(with (3 <C 1) simultaneously implies a smallness of the "self-field" effect and vice versa. 
And since artificially prepared Josephson networks allow for more flexibility in vary- 
ing the experimentally-controlled parameters, it is always possible to keep the both 
above effects down by appropriately tuning the ratio d/\j. Typically, in this kind of 
experiments S = 0.01 — O.lfim 2 and d <C \j. 

In what follows, we are interested in the magnetic field induced behavior of the 
electric polarization in a 3D JJA at zero temperature. Recall that a conventional 
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(zero-field) pair polarization operator within the model under discussion reads [30] 



N 

P = J2wi, (16) 
1=1 

where qi = — 2en,i with rii the pair number operator, and r\ is the coordinate of the 
center of the grain. 

In view of Eqs.(l)-(6), and taking into account a usual " phase- number" commuta- 
tion relation, [(pi,nj] = i5 it j, the evolution of the pair polarization operator is deter- 
mined via the equation of motion 

f = ^ \P,n] = |E Jsm<f> hJ (H)r h3 (17) 

ij 

Resolving the above equation, we arrive at the following net value of the magnetic-field 
induced polarization (per grain) 

P(H) = -< p{t) > = —j dtj dt>Y,< BinMmj >, (18) 

o o v 

where < ... > denotes a configurational averaging over the grain positions, while the 
bar means a temporal averaging (with a characteristic time r). To consider a field- 
induced polarization only, we assume that in a zero magnetic field, P = 0, implying 

<Mo) = o. 

Since, as usual, the main qualitative results presented here do not depend on the 
particular choice of the probability distribution function, for a change, the following 
law f(f,R) = f r (r)f R (R) with f r (r) = ^^e-^^l^ and f R (R) = 5(X - 
d)5(Y — d)5(Z — d) will be assumed. As a result, we observe that the magnetic field 
H = (0, 0, H z ) (applied along the z-axis) will induce a non-vanishing longitudinal (along 
a;-axis) electric polarization 

P X (H Z )=P G(H Z /H ), (19) 

with 

G(z) = ze-" 2 (20) 

Here Po = edrJ/h, H = §o/S with S = ird 2 being an average projected area of a 
single junction, and z = H z /H . For small applied fields (z <C 1), the induced polar- 
ization P X (H Z ) m aoH z , where cto = erdJ/hH is the so-called linear magnetoelectric 
coefficient [35]. However, as we mentioned in the very beginning, to correctly describe 
any induced effects in very small external fields, both Josephson junction pinning and 
"self-field" effects have to be taken into account. At the same time, in view of Eq.(16) 
the induced polarization is related to the corresponding change of the effective charge 
5Q = (l/d) J2i < <li x i > in applied magnetic field as follows 

5Q(H Z ) = ^j^- = Q G(H Z /H ), (21) 
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Figure 2: The behavior of the induced effective charge SQ/Q in applied magnetic field 
H z /H , according to Eq.(21). 

where Q = erJ/h. 

It is of interest also to consider the related field behavior of the effective flux ca- 
pacitance 5C = Td5Q(H z )/d§ which in view of Eq.(21) reads 

SC(H Z ) = C (l - 2%^ e -^/H^ (22) 

where <£> = SH Z , and Co = tQ /§ . Figure 2 shows the behavior of the induced 
effective charge 5Q/Q in applied magnetic field H z /H . As is seen, at H z /H 0.8 
the effective charge reaches its maximum (while the capacitance changes its sign at this 
field), suggesting a significant redistribution of the junction charge balance in a model 
system under the influence of an applied magnetic field, near the Josephson critical 
field H . Note that a somewhat similar behavior of the magnetic field induced charge 
(and related capacitance) has been observed in 2D electron systems [36]. 

Taking r = ICr 10 s for the Josephson relaxation time (which is related to the 
Josephson plasma frequency, oo p ~ r _1 , known to be the characteristic frequency 
of the system for Ec ^ J regime, with Ec being the Coulomb grain's charge en- 
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ergy; typically u p = 10 9 — 10 11 Hz) and J/ks = 90K for a zero-temperature Joseph- 
son energy in YBCO ceramics, we arrive at the following estimates of the effec- 
tive charge Qo ~ 10~ 16 C, flux capacitance Cq ~ 10 _11 -F, the equivalent current 
-^o ~ Qo/t ~ 10~ 6 A, and voltage V ~ Qq/Cq lO -5 ^. We note that the above 
set of estimates fall into the range of parameters used in typical experiments to study 
the charging effects both in single JJs (with the working frequency from RF range of 
uj ~ 10GHz used to stimulate the system [13]) and JJAs [37] suggesting thus quite an 
optimistic possibility to observe the above-predicted field induced effects experimen- 
tally, using a specially prepared system of arrays of superconducting grains. 

4. STRESS INDUCED PARAMAGNETIC MEISSNER EFFECT 

The possibility to observe tangible piezoeffects in mechanically loaded grain bound- 
ary Josephson junctions (GBJJs) is based on the fact that under plastic deformation, 
grain boundaries (GBs) (which are the natural sources of weak links in HTCS), move 
rather rapidly via the movement of the grain boundary dislocations (GBDs) comprising 
these GBs [38-40]. Using the above evidence, in Ref.8 a piezophase response of a single 
GBJJ (created by GBDs strain field acting as an insulating barrier of thickness I and 
height U in a SIS-type junction with the Josephson energy J oc e~ 1 ^) to an applied 
stress was considered. To understand how piezoeffects manifest themselves through 
GBJJs, let us invoke an analogy with the so-called thermophase effect [6,7] (a quantum 
mechanical alternative for the conventional thermoelectric effect) in JJs. In essence, the 
thermophase effect assumes a direct coupling between an applied temperature drop AT 
and the resulting phase difference A<fi through a J J. When a rather small temperature 
gradient is applied to a JJ, an entropy-carrying normal current /„ = L n AT (where L n 
is the thermoelectric coefficient) is generated through such a junction. To satisfy the 
constraint dictated by the Meissner effect, the resulting supercurrent I s = / c sin[A0] 
(with I c = 2eJ/h being the Josephson critical current) develops a phase difference 
through a weak link. The normal current is locally canceled by a counterflow of super- 
current, so that the total current through the junction I — I n + I s — 0. As a result, 
supercurrent / c sin[A0] = — I n = —L n AT generates a nonzero phase difference leading 
to the linear thermophase effect [6,7] A0 ~ —L tp AT with L tp = L n /I C (T). 

By analogy, we can introduce a piezophase effect (as a quantum alternative for the 
conventional piezoelectric effect) through a J J [8,16]. Indeed, a linear conventional 
piezoelectric effect relates induced polarization P n to an applied strain e as [35] P n = 
d n e, where d n is the piezoelectric coefficient. The corresponding normal piezocurrent 
density is j n = dP n /dt = d n e where e(<r) is a rate of plastic deformation which depends 
on the number of GBDs of density p and a mean dislocation rate Vd as follows [41] 
e(cr) = bpvd((j) (where b is the absolute value of the appropriate Burgers vector). In 
turn, Vd(o~) ~ vo(a/a m ). To meet the requirements imposed by the Meissner effect, in 
response to the induced normal piezocurrent, the corresponding Josephson supercurrent 
of density j s = dP s /dt = j c sin[A0] should emerge within the contact. Here P s = 
—2enb is the Cooper pair's induced polarization with n the pair number density, and 
j c = 2ebJ/hV is the critical current density. The neutrality conditions (j n +js = and 
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P n + P s = const) will lead then to the linear piezophase effect A0 ~ —d pp e(o) (with 
d pp = d n /j c being the piezophase coefficient) and the concomitant change of the pair 
number density under an applied strain, viz., An(e) = d pn e with d pn = d n /2eb. 

To adequately describe magnetic properties of a granular superconductor, we again 
employ a model of random three-dimensional (3D) overdamped Josephson junction ar- 
ray which is based on the familiar tunneling Hamiltonian given by Eq.(l) (see Section 
2). According to the above-discussed scenario, under mechanical loading the supercon- 
ducting phase difference will acquire an additional contribution 5<pi t j(o) = —Bo ■ r*jj, 
where B = d n eo/o m j c b with eo = bpvo being the maximum deformation rate and the 
other parameters defined earlier. If, in addition to the external loading, the network of 
superconducting grains is under the influence of an applied frustrating magnetic field 

— * 

H, the total phase difference through the contact reads 

<M#, a) = <Pl + ^(f itj A R hj ) ■ H — Bo ■ r h3 (23) 

Once again, to neglect the influence of the self-field effects in a real material, the 
corresponding Josephson penetration length Aj must be much larger than the junction 
(or grain) size. Likewise, to ensure the uniformity of the applied stress o, we also assume 
that d <C A CT , where A CT is a characteristic length over which o is kept homogeneous. 

When the Josephson supercurrent Jf ■ = I c sin faj circulates around a set of grains, 
it induces a random magnetic moment p s of the Josephson network which results in 
the stress induced net magnetization (Cf. Eq.(8)) 

oo 

M S (H, a) =< p s >= J dnjdRijfinj, R i:j )p s (24) 

o 

To capture the very essence of the superconducting piezomagnetic effect, in what follows 
we assume for simplicity that an unloaded sample does not possess any spontaneous 
magnetization at zero magnetic field (that is M s (0, 0) = 0) and that its Meissner 
response to a small applied field H is purely diamagnetic (that is M s (H,0) ~ —H). 
According to Eq.(23), this condition implies 0° ■ = 2irm for the initial phase difference 
with m = 0,±1,±2, ... 

Taking the applied stress along the x-axis, o = (o, 0,0), normally to the applied 
magnetic field H = (0,0, H), and assuming an exponential distribution law for the 
distance between grains, f r (f) = f(x)f(y)f(z) with f(xj) = (l/d)e~ x ^ d (see Section 2 
for detailes), we get finally 

M,(.,^-M 0W [1 ^)/^ r ( 25 > 

for the induced transverse magnetization, where H tot (H,o) = H — H*(o~) is the total 
magnetic field with H*(o) = [o / o (o)]H being a stress-induced contribution. Here, 
M (o) = I c (o)SN/V with S = ndD being a projected area around the JJ, H = &o/S, 
and er (er) = o m [j c (o) / jd\ip/d) with j d = d n e and e = bpv being the maximum 
values of the dislocation current density and the plastic deformation rate, respectively. 
According to the recent experiments [43], the tunneling dominated critical current 
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Figure 3: The reduced magnetization M s /M c as a function of the reduced applied 
magnetic field H/H , according to Eq.(25) for different values of reduced applied stress: 
a/a c = (solid line), a/a c = 0.01 (dotted line), and cr/a c = 0.05 (dashed line). 

I c (and its density j c ) in HTCS ceramics was found to exponentially increase under 
compressive stress, viz. / c (o") = i c (0)e /3cr with f3 ~ \jo m . Specifically, the critical 
current at a — 9kbar was found to be three times higher its value at a = l.Bkbar, 
clearly indicating a weak-links-mediated origin of the phenomenon (in the best defect- 
free thin films this ratio is controlled by the stress induced modifications of the carrier 
number density and practically never exceeds a few percents [42]). Strictly speaking, 
the critical current will also change (decrease) with applied magnetic field. However, 
for the fields under discussion (see below) in the first approximation this effect can 
be neglected. In view of Eq.(25), dependence of I c on a will lead to a rather strong 
piezomagnetic effects. Indeed, Fig. 3 shows changes of the initial stress-free diamagnetic 
magnetization M s /M c (solid line) under an applied stress o jo c . Here M c = M (0) and 
o c = o" (0) (see below for estimates). As we see, already relatively small values of 
an applied stress render a low field Meissner phase strongly paramagnetic (dotted and 
dashed lines) simultaneously increasing the maximum of the magnetization and shifting 
it towards higher magnetic fields. According to Eq.(25), the initially diamagnetic 
Meissner effect turns paramagnetic as soon as the piezomagnetic contribution H*(a) 
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exceeds an applied magnetic field H. To see whether this can actually happen in a 
real material, let us estimate a magnitude of the piezo magnetic field H*. Typically, for 
HTCS ceramics S ~ 10/im 2 , leading to H$ ~ 1G. To estimate the needed value of the 
dislocation current density jd, we turn to the available experimental data. According 
to Ref.39, a rather strong polarization under compressive pressure a/a rn ~ 0.1 was 
observed in YBCO ceramic samples at T = 77 K yielding d n = 10 2 C/m 2 for the 
piezoelectric coefficient. Usually, for GBJJs eo — 10~ 2 s _1 , and b ~ lOnm leading to 
3d — d n e ~ lA/m 2 for the maximum dislocation current density. Using the typical 
values of the critical current density j c (cr) = 10 4 A/m 2 (found [43] for a/a m ~ 0.1) 
and grain size d ~ 1/im, we arrive at the following estimate of the piezomagnetic field 
H* ~ 10~ 2 H . Thus, the predicted stress induced paramagnetic Meissner effect should 
be observable for applied magnetic fields H ~ 10~ 2 H ~ 0.01G which correspond to 
the region where the original PME was first registered [1-3]. In turn, the piezoelectric 
coefficient d n is related to a charge Q in the GBJJ as [44] d n = (Q/S)(d/b) 2 . Given the 
above-obtained estimates, we get Q ~ 10~ 13 C for an effective charge accumulated by 
the GBs. Notice that the above values of the aplied stress a and the resulting effective 
charge Q correspond (via the so-called electroplastic effect [44]) to an equivalent applied 
electric field E = b 2 a/Q ~ 10 7 U/m at which rather pronounced electric-field induced 
effects in HTCS have been recently observed [17-21]. 

Besides, according to Ref.43 the Josephson projected area S slightly decreases un- 
der pressure thus leading to some increase of the characteristic field H = &o/S. In 
view of Eq.(25), it means that a smaller compression stress is needed to actually reverse 
the sign of the induced magnetization M s . Furthermore, if an unloaded granular super- 
conductor already exhibits the PME, due to the orbital currents induced spontaneous 
magnetization resulting from an initial phase difference 0^ = 2nr in Eq.(23) with 
fractional r (in particular, r = 1/2 corresponds to the so-called [2,3,21] 7r-type state), 
then according to our scenario this effect will either be further enhanced by applying 
a compression (with a > 0) or will disappear under a strong enough extension (with 
a < 0). Given the markedly different mechanisms and scales of stress induced changes 
in defect-free thin films [42] and weak-links-ridden ceramics [43], it should be possible 
to experimentally register the suggested here piezophase effects. 
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